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OPTIMAL STRATEGIES OF INVESTMENT IN A LINEAR 
STOCHASTIC MODEL OF MARKET 

G. S. KAMBARBAEVA AND O. S. ROZANOVA 


Abstract. We study the continuous time portfolio optimization model on the 
market where the mean returns of individual securities or asset categories are 
linearly dependent on underlying economic factors. We introduce the functional 
featuring the expected earnings yield of portfolio minus a penalty term pro¬ 
portional with a coefficient 7 to the variance when we keep the value of the 
factor levels fixed. The coefficient 7 plays the role of a risk-aversion parame¬ 
ter. We find the optimal trading positions that can be obtained as the solution 
to a maximization problem for at any moment of time. The single-factor 
case is analyzed in more details. We present a simple asset allocation example 
featuring an interest rate which affects a stock index and also serves as a sec¬ 
ond investment opportunity. We consider two possibilities: the interest rate for 
the bank account is governed by Vasicek-type and Cox-Ingersoll-Ross dynam¬ 
ics, respectively. Then we compare our results with the theory of Bielecki and 
Pliska where the authors employ the methods of the risk-sensitive control the¬ 
ory thereby using an infinite horizon objective featuring the long run expected 
growth rate, the asymptotic variance, and a risk-aversion parameter similar to 
7 - 


1. Introduction 

The art of making decisions about investment mix in order to meet specified 
investment goals for the benefit of the investors and balancing risk against perfor¬ 
mance is called the portfolio management (portfolio is a collection of investments all 
owned by the same individual or organization). 

A modern study of portfolio selection begins in works by Markowitz [55] , [55]. He 
showed how to formulate the problem of minimizing a portfolio’s variance subject 
to the constraint that its expected return equals a prescribed level as a quadratic 
program. Such an optimal portfolio is said to be variance minimizing and if it 
also achieves the maximum expected return among all portfolios having the same 
variance of return then it is said to be efficient m- 

There has been considerable research involving stochastic processes models of 
assets taking into account optimal investment decisions. Several researchers (e.g. 
Merton (31], Karatzas |25|) used stochastic control theory to develop continuous 
time portfolio management model where the assets are modeled bas stochastic pro¬ 
cesses but financial and economic factors are ignored. At the same time, a number 
of empirical studies (e.g. [31], [33], [5D]) provided evidences that macroeconomic fac¬ 
tors such as unemployment rate, inflation rate, dividend yield, change of industrial 
production, an interest rate, etc, influence on the stock return. 

Lucas m introduced a discrete time model including stochastic process models 
as factors. Brennan, Schwartz and Lagnado [7] considered the factors as diffusion 
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processes and assets as correlated Brownian motions, with the drift and diffusion 
coefficients for the asset pricesses taken to be deterministic functions of the factor 
levels. Their objective was to maximize expected utility of wealth at a terminal date. 
A key limitation of the Brennan-Schwartz-Lagnado approach is that one is unlikely 
to obtain tractable formulas for the optimal strategies. 

In the past decades, the applications of risk-sensitive control to asset management 
is very popular. The risk-sensitive control differs from traditional stochastic control 
in that it explicitly models the risk-aversion of the decision maker as an integral part 
of the control framework, rather than importing it in the problem via an externally 
defined utility function [41) . Risk-sensitive control was first applied to solve finan¬ 
cial problems by Lefebvre and Montulet |26j in a corporate finance context and by 
Fleming m in a portfolio selection context. 

Bielecki and Pliska [3] were the first to apply the continuous time risk-sensitive 
control as a practical tool that could be used to solve ’’real world” portfolio selection 
problems. In the series of works by T.Bielecki, S.Pliska et al.(|4], [5], etc) was devoted 
to a infinite-horizon continuous-time risk sensitive portfolio optimization problem. 
The authors considered a model of market analogous to the Brennan-Schwartz- 
Lagnado one [7] where the mean returns of individual securities are explicitly affected 
by underlying economic factors such as dividend yields, a firm’s return on equity, an 
interest rate, and unemployment rate. The factors are random processes, and the 
drift coefficients for the securities are linear functions of these factors. The main 
result of the theory by Bielecki and Pliska is a construction of admissible trading 
strategies, which have a simple characterization in terms of the factor levels. The 
results are illustrated on a simple but important example of two asset allocation, 
having independent interest to financial economists. Here one of assets is a bank 
account and the unique factor is a Vasicek-type interest rate. The Vasicek model of 
the interest rate is linear and this gives a possibility to obtain an explicit formulas 
for the optimal strategy. 

The strategy proposed in the works of Bielecki and Pliska refers to the strategic 
asset allocation. According the primary goal of a strategic asset allocation is to 
create an asset mix that will provide the optimal balance between expected risk and 
return for a long-term investment horizon. 

In the present paper we propose an alternative method of capital allocation in 
which an investor takes a more active approach that tries to position a portfolio 
into those assets, sectors, or individual stocks that show the most potential for 
gains. The model refers to the tactical asset allocation (e.g. [35]). Our strategy also 
has a simple characterization in terms of the factor levels, but it is more flexible 
comparing with the Bielecki-Pliska model and can be actualized within all the time 
of investment. Moreover, we get explicit formulae not only for the lineal model 
of factor (in particular, for the Vasicek model of the interest rate), but for more 
complicated one such that the Cox-Ingersoll-Ross model. This paper summarizes 
and extends the results of [21], [22], [23] . 

This paper is organized as follows. In Sec[^ we describe the model of linear 
market where we are going to consider the allocation of capital. In Sec[^we give an 
outline of the Bielecki and Pliska theory and write their explicit optimal strategy. 
Sec]^ contains auxiliary results: an algorithm of finding the conditional expectation 
and conditional variance for a couple of stochastic differential equations. To find 
these values we have to solve a Fokker-Planck equation for the joint probability 
density of two respective random values. We show that there exist two approaches 
to solving this problem. The first approach uses an ansatz for the solution, such 
that the problem is reduced to solving a nonlinear system of ordinary differential 
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equations. The second approach refers to the application of the Fourier analysis. 
In Secj^ we formulate the problem of the fixed time portfolio selection and give a 
general algorithm for its solution. 

In Secj^we consider the case of linear interest rate (the Vasicek model). First we 
solve the fixed time optimization problem for the the example of portfolio consisting 
of two assets mentioned before. Then we find asymptotics of the proportions of 
capital invested in securities as the time tends to infinity for different initial distri¬ 
bution of factor. We dwell on the case of two and three risky assets and discuss the 
influence of parameters of the model on the strategy. At last we compare our asset 
allocation with its analogous for the long-run Bielecki and Pliska strategy. 

Secj^ deals with the Cox-Ingersoll-Ross model of the interest rate, we consider 
again the case of portfolio consisting of two assets. 

In Secj^we compare the fixed time optimal strategies for a portfolio consisting of 
two assets for both model of interest rate and conclude that the Cox-Ingersoll-Ross 
model is in some sense preferable. 

The formulae appearing in this work are sometime very cumbersome and we have 
no possibility to write them out. To obtain them we used the computer algebra 
system MAPLE. 


2. A STOCHASTIC MARKET MODEL 

We study a portfolio optimization problem in the frame of the market model of 
m >2 assets and n > I factors used by T.Bielecki and S.Pliska (e.g. mM)- Below 
we describe this model. 

Let (n, {J^t}t>o, .A, P) be the underlying probability space. Denoting by Si, i = 
the price of the i-th security and by Xj, j = the level of the j-th 

factor at time t, we consider the following market model for the dynamics of the 
security prices and factors: 


^ — (Ai -t- ^ aipXp{t))dt 

' p—1 

Si{0) = Si > 0,i = 1, 


m-\-n 

+ ^ cTikdWkit), 


( 2 . 1 ) 


m+n 


dXj{t) = {Bj + ^j3jpXp{t))dt+ ^ \jkdWk{t), 

p=i 




( 2 . 2 ) 


where W (t) is a 7^™+’^ - valued standard Brownian motion process with components 
Wk{t)-, X{t) is the 7^" - valued factor process with components Xj{t)-, the market 
parameters A := [A*], B := [Bj\, a := [aip], j3 := [Pjp], E := [aik], A := [Ajfc] 
matrices of appropriate dimensions. According to [24] (chapter 5) a unique, strong 
solution exists for (2.1), (2.21, and the processes Si{t) are positive with probability 1. 

Let Qt ■= (^{{S{s),X{s)),0 < s < t), where S{t) = {Si{t),Sm{t)) is the se¬ 
curity price process. Let h{t) = denote an 7?.™ valued investment 

process or strategy whose component hi{t) represents the proportion of capital that 
is invested in security i at time t. We define the admissible investment strategy 
according to |4]. 
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Definition 2.1. An investment process h{t) is admissible if the following conditions 
are satisfiecQ 

m 

(*) = 

2 — 1 

(ii) h(t) is measurable, C/t—adapted; 

(Hi) P[ / k’"{s)h{s)ds < oo] = 1 for all finite t > 0. 

Jo 


The class of admissible investment strategies will be denoted by H. 

Let h{t) be an admissible investment process. Then there exists a unique, strong 
and almost surely positive solution V(t) to the following equation: 


m 

dV{t) =Y,K{t)V{t) 
v{0) = v>0. 


n m+n 

{Ai + ^ aipXp{t))dt + ^ (TikdWk{t) , 
p=i fe=i 


(2.3) 


The process V{t) represents the investor’s capital at time t where hi{t) represents 
the proportion of capital that is invested in security i. 


Remark 2.1. In [T3] the model of the linear market was extended to the case of 
asset prices represented by SDEs driven by Brownian motion and a Poisson random 
measure, with drifts that are functions of an auxiliary diffusion factor process. 


3. The optimal investment strategy by Bielecki and Pliska [1] 


A new kind of portfolio optimization model to the type of asset allocation prob¬ 
lems (2.1) considering a portfolio of m > 2 assets affected by n > 1 financial and 
economic factors was introduced by Bielecki and Pliska. Namely, they considered 
the following functional 

Qeit) 


Jg := liminf ■ 


t 


Qs{t) ■■= -^lnE( 


,(-e/2)lnV(t) 


> - 2 , 


A Taylor expansion of Qg around 9 — 0 yields 


Qg{t) = E(lnP(t)) - -Var(lnP(t)) + 0(0^), 


(3.1) 


hence Jg can be interpreted as the long-run expected growth rate minus a penalty 
term, with an error that is proportional to . The penalty term is also proportional 
to 0, so 9 was interpreted as a risk sensitivity parameter or risk aversion parameter, 
with 0 > 0 and 0 < 0 corresponding to risk averse and risk seeking investors, 
respectively and 0 = 0 is the risk null case. 

Bielecki and Pliska [4] proposed to solve the following family of risk sensitive 
optimal investment problems, labeled as {Pg)'. 

for 9 G (0,oo) maximize the risk sensitized expected growth rate 

Je(u,x;h(-)) = liminf =u,A(0) =a:] 

t—>-oo 0 


over the class of all admissible investment processes h(-), subject to 
definition 1.3, where V{t),X{t) obey equations \2.2\ . 


^(•)^ stands for a transposition operator 

^E(-) and Var{-) are expectation and variance in the probability space P), 

respectively 
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Authors noticed that Jg has a large-deviations-type functional for the capital 
process V (t) . Maximizing Jg for 9 > 0 protects an investor interested in maximizing 
the expected growth rate of the capital against large deviations of the actually 
realized rate from the expectations. 

Remark 3.1. In case of 0 < 0 the problem (Pg) can be solved similar to the case of 
6 > 0. The case 0 = 0 is considered separately as a limiting case 0 —)■ 0. 


An algorithm to find optimal investment strategy Hg together with corresponding 
maximum value of Jg labeled as p{9) was proposed in [3]. In order to present 
the main results pertaining to these investment problems, authors introduced the 
following notation for 0 > 0 and x G i?”: 


Kg(x) 


inf 


^ + 1^ {A + ax) 


(3.2) 


where {A + ax) denotes the vector with components (A + ax)i = [Ai + E (XipXp^. 
Also they made the following assumptions: 


Assumption 3.1. x = 

Assumption 3.2. lim Kg{x) = —oo for 0 > 0. Here | 
IkIKoo 

Assumption 3.3. The matrix AA"^ is positive definite. 
Assumption 3.4. The matrix EA^ is zero. 


is the norm in K". 


Remark 3.2. (i) Note that if EE^ is positive definite, then assumption 3.2 is implied 
by assumption |3.1| 

(ii) Assumptions 3.1|3.4 are sufficient for the results below, but Assumption 3.2 
is not necessary (see [5], Sec. 4). 

Following two theorems contain key results concerning the solution of the problem 

(Po)- 

Theorem 3.1. [3] Assume 3. IE ^ For a fixed value of 9 > 0 let Hg{x) denote a 
minimizing selector in (3.2), that is 


Kg{x) := 


^ ^ + Ij Hg{x)'^Y.Y.^Hg{x) - Hg{x)^{A + ax) 


Then the investment process hg is optimal for problem (Pg), where Vt > 0 

he{t) = HgiXit)). 


(3.3) 


Theorem 3.2. Let us assume \3.1^3.3\ and consider problem (Pg) for a fixed value 
of 9 > 0. Let hg(f) satisfy theorem 3.1. Then 

(1) For all V > 0 x G i?" we have 

Jg{v,x-,hg{-)) = lim =u,A(0) =:r] 

i->oo \ 9 J 

■■= p{0)- 
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(2) The constant p{9) is the unique non-negative constant which is a part of the 
solution {p{9),v{x;9)) to the following equation: 


P 


{B + j3x)^ gradxvix) 





dv{x) dv(x) 
dxi dxj 


n+m 






1 a2 / \ 

+ 2 ^ ^ KkXjk-Keix), 

iJ — 1 k—1 

v{x) G C^(ii”), lim v(x) = oo, p = const. 
Ikll-s-oo 


(3.4) 


It remains to consider the case corresponding to 0 = 0. This is the classical 
problem of maximizing the portfolio’s expected growth rate, that is, the growth rate 
under the log-utility function (see e.g. Karatzas [IS]). This problem was labeled 
(Po) and formulated as follows: 

maximize the functional 

Jq{v,x-, h{-)) = liminff^ lnE[lny(t)|T/(0) = u,X(0) = x] 


over the class of all admissible investment processes h{-) subject to 
definition 1.3, where V{t),X(t) described by equations (2.3), (2.2). 

It turns out that to solve (Pq) h is necessary to make three additional assumptions: 


Assumption 3.5. For each 9 > 0 the function Kg{x) defined in (3.2) is of the quadratic 
form 

Kg(x) = ^x'^Ki{9)x -f K2{9)x -k K3{9), 


where Ki{9), K2{9) Kz{9) are functions of appropriate dimensions depending only 
on 9. 


Assumption 3.6. For each 0 > 0 the matrix Ki{9) is symmetric and negative definite. 


Assumption 3.7. The nx n matrix /3 with components Pjp in (2.2) is stable. 

is satisfied if, for example, the matrix is non- 


Remark 3.3. (i) Assumption 
singular and y = i?". 


3.5 


(ii) According to |3| Assumption 3.5 implies lim Ki{9) = Ki{Q) for i = 1, 2,3. 


e-fo 


In order to establish relationships between the risk-neutral problem (Pq) and the 
risk-sensitive problem (Pg), 0 > 0 we consider the following equation: 


p(0) = (P -I- Px)'^gradxVo{x) -|- 


1 u2 / ^ n+m 

1 d^voix) , 

2^ XikXjk Ao(a 


vo{x) e C\R^), 


2 dxidxj 
1,3 = 1 •> 

lim V{ 3 {x) = 00 , 


k=l 

p(0) = const. 


(3.5) 


The following two results are true: 

Theorem 3.3. |3| Assume |5'.5]|5'. 7] Then the optimal strategy for (Pg) is as in 
Theorem 3.1 with 9 — 0, and the optimal objective value p{0) = p{9) in Theorem 3.2 
with 9 = 0. Moreover, the optimal objective values p(9),9 > 0, being the solution of 
(3.4), converge to the optimal objective value p{0) as 9 ^ 0. 


The next result characterizes the portfolio’s expected growth rate under the op¬ 
timal investment strategy for the risk aversion level 0 > 0. We denote this growth 
rate by pg, which is to be distinguished from the optimal objective value p{9), as in 
Theorem 13.11 
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Theorem 3.4. [1] Assume 3.3^3 . 7 Fix 0 > 0, let Hg{x) be as in theorem \3.1\ and 
suppose that Fig (a;) is an ajfine function and that 

lim — Hg{x)^{A + aa;)] = —oo. 

||x||—>00 2 


Consider the equation 

. .rp 1 d'^vgo{x)'^^ 

pg = {B + (3xY grad:^vg.Q{x) + a ^ ^ ^ 


-[\Hg{x) 

Vg,o{x) e 


2 ^ dxidxj 
1,3 = 1 ■> 


k=l 

\T( 


Hg{x) — Hg{x)^ {A + ax)], 
lim vg^o(x) = cx), pg = const. 


(3.6) 


Then there exists a solution {pg,vgfi) the preceding equation, the constant pg is 
unique, and we have 

Jo{v,x;hg{-)) = pg 

for all {v,x) S ( 0 , 00 ) x i?", where hg(-) is defined as in (3.3). 

The main result of Bielecki and Pliska is that the optimal investment strategy 
problem is converted to the problem of solution to PDE (3.4) . Authors solve the 
problem explicitly for the classical example of portfolio consisting of two assets, 
where one of them is a bank account, and a linear interest rate as a factor. 

Namely, they consider a single risky asset, say a stock index, that is governed by 
a stochastic differential equation 
dSi{t) 


Slit) 


= {Ai + aiRit))dt + aidWift), S'i(O) = s > 0, 


where the spot interest rate R{t) is satisfies the classical Vasicek dynamics: 

dR{t) = {B + l3Rit))dt + XdW 2 {t), i?(0) = r. 

Here Ai, ai, H,/3, CTi, A are fixed, scalar parameters to be estimated, while Wi,W 2 
are two independent Brownian motions. Hereafter we assume H > 0,/3 < 0. in all 
that follows. 

The investor can take a long or short position in the stock index as well as borrow 
or lend money, with continuous compounding, at the prevailing interest rate. It 
is therefore convenient to follow the common approach and introduce the ’’bank 
account” process S' 2 , where 

dS2{t) 


S2it) 


= Rit)dt. 


Thus S' 2 (t) represents the time t value of a savings account when 52(0) = 1 dollar 
is deposited at the zero time. 

With only two assets it is convenient to describe the investor’s trading strategy 
in terms of the scalar valued function Hg which is interpreted as the proportion of 
capital invested in the stock index, leaving the proportion 1 — Hg invested in the 
bank account. 

This enables us to formulate the investor’s problem as [Pg) in the market model 
( plj ), for there are m = 2, n = 1, and we can set (/ii,/i 2 ) = [Hg,! — Hg), 
X{t) = Rft), H = H, /3 = ^, A = (0,0, A)^, A = (Ai, 0)^, a = (ai, 1)^, 

cTi 0 0 

^000 

According to Theorem |3.1| 


E = 


Kg(R) = inf [(l/2)(6»/2 + l)(h, 1 - h)EE^(/i, 1 - h)^ - [h, 1 - h){A + aR)], 
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in which case 


He := He{R) = 


Ai + (cki — 1)R 


(1 + 1 ) 




Ke{R) = -R- 


(^1 + (gi ~ 

(0 + 2)c 


Theorem 3.2 implies that p(0) is the part fo solution {p,v) of the equation 
p = ^X^v"iR) + {B + m)v\R) - \\\v'{R)Y - Ke{R), 
and, according to [5], equals 


\ 2/3 ]\t2 

p{6) = A^iVi + BN2 - 


Al 


(0 + 2), 


r2 ’ 


where 


Q I I rp” ^ 2Ai{a.i — l) o p A,r 

,, p + yp + ie+ 2 )ai 1- ie+ 2 )ai + 2SiVi 

-‘''1 — ”—iV2 — ■- 


A20 


02 I e\^{ai-iy 

P (B+ 2 )aj 


Theorem |3.4| results that the solution of equation 
1 

2 ' 


p = lX'^v''{R) + {B + m)v'{R)- 


- ^(HeiR), 1 - He{R))X:XA{He{R), 1 - He{R)V- 
-{HgiR),l- Hg{R)){A + aR)], 


gives 


B 


P0 = 


2(0 + 1 ) 


/3 (0 + 2)V2 L 




2/3 


(3.7) 


(3.8) 


(3.9) 


We note that Bielecki and Pliska introduced an optimal investment strategy to 
maximize portfolio return to an infinite time horizon. In this paper we introduce 
another strategy which can be used by investor to manage portfolio and maximize 
return at any fixed time moment. 


4. Conditional expectation and variance for a couple of SDEs: two 

APPROACHES TO THE SOLUTION 
Let us consider a system of stochastic differential equations 
dF = A{t,F,X)dt + a{t,F,X)dWi, 

dX = B{t,F,X)dt +X{t,F,X)dW2, (4.1) 

F(0) = /, X(0) =x,t>0,feR,xeR, 

where W = (Wi, 1+2) ia a two-dimensional Brownian motion with independent com¬ 
ponents, A, B, a, X are given functions. 

The joint distribution density P{t,f,x) of stochastic variables F and X is de¬ 
scribed by the Fokker-Plank equation (e.g. [35], |33] 1 


dP{t,f,x) 

dA{t,F,X)P{t,f,x) ld^a\t,F,X)Pit,f,x) 


dt 

df 2 dp 

dB{t,F,X)PitJ,x) ldPpt,F,X)Pit,f,x) 

dx 2 dx"^ 

(4.2) 

with initial data 

P(0,/,x) = Poif,x), 

(4.3) 


determined by initial distributions of F and X. 
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Provided P{t,f,x) is known, the conditional expectation of F with given value 
of X at the moment t can be found by the following formula (see, e.g. m, 0) 


fit,x) :=-E{F\X = x) = 




(4.4) 


If we set P[j{f,x) = 6{f — fo)g{x), where /o S M and g{x) is an arbitrary function, 
Jg^g{x)dx = 1, then /(0,x) = /q. Some characteristics of (4.41 are studied in [T], [5]. 

The conditional variance of a stochastic variable F with given value of X at the 
moment t is defined as follows: 


~i{t,x) := Var(F|X = x) = 


iRpPjt, f,x)df 

PPitJ^pdf 




(4.5) 


The fundamental solution of equation ( |4.2[ ) can be found by means of the Riccati 
matrix equations ga, El. For some simple but important for application choice of 
initial data problem ( |4.2 1, ( |4.3[ ) can be solved in terms of elementary functions. 

Moreover, sometimes the Fourier transform of P{t, /, x) can be found easier than 
this function itself. Further we are providing two approaches to the problem. 


4.1. Approach 1: a reduction to the system of ODEs. Let us assume that 
the variables F and X obey the following system of stochastic differential equations: 


dF{t) = (A + aiX{t) + a 2 F{t))dt + aidWi{t) + a 2 dW 2 {t), 
dX{t) = {B + (}\X{t) + j32F{t))dt + XidWi{t) + X2dW2{t), 
F{0)=f,X{0)=x,t>0,f,xeR, 


(4.6) 


where W{t) = {Wi{t),W 2 {t)) is a two-dimensional Brownian motion; A, B, ai,/3i,ai, Xi are 
known smooth functions oi t, i = 1, 2. 

Then joint distribution density P{t,f,x) of F and X solves the Fokker-Planck 
equation 


dP{t,f,x) dP{tJ,x) 

- + ai{t)x) - — - 

-02 {t) (P{t, f,x) + f 


-/3i(i)(P(t,f,x) 


df 

dP{t,f,x) 
dx 


-{B{t)+^2it)f) 
1 


dP{t,f,x) 


) + + Pit)) 


dx 


(4.7) 


-f 


+ (o’i(t)Ai(t) -I- cr2{t)X2{t)) -^ + p{t)) 


dp 


dx^ 


subject to initial data 


P{0,f,x) = Po{f,x) = S{f - fo)g{x). 


(4.8) 


We perform the Fourier transform of P{t,f,x) in / of (4.7) and (4.8) and get: 
dP{t, /i, x) 


dt 


= —{A{t) + ai{t)x)gP(t, g, x)i— 


-a 2 {t) P{t,g,x) - g 


dPit, g, x) 
dfi 


- {B{t) + di{t)x) 


dP{t, fi, x) 
dx 


-l3i{t)P{t,g,x) - I32it) ^ ^dxd^^'^ ^~ 

+ Pit))pPit,g-,x) + ^(A?(t)-f A^(t))^-^|^- 
-|-(cri(t)Ai(t) -I- a 2 (t)X 2 it))g, ^^^^g^’^^ i, 


(4.9) 
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P(0,/i,a;) = 


1 






9{x)- 


To obtain explicit formulae we restrict ourselves to the case 

^- 

9{x) = 


(4.10) 


(4.11) 


where ccq € K, s S M+. Here a;o is the mean of the variable X at the initial time, 
is the variance. Thereafter we seek the solution of the problem (4.9), (4.10) using 
the following ansatz: 


P{t,n,x) = 


D7i(^)+72(i)M+73(i)a:+74Ci)M +75(i)a:M+76(i)a^ 

svTtt 


(4.12) 


We substitute \A.12\ into ( |4.9[ ) and (4.10). Equating coefficients at the powers of 
p, and X gives the system for = 1,..., 6: 


dt 

dt 


dlsjt) 

dt 

dli{t) 

dt 

dibit) 

dt 


dl&jt) 

dt 


imilit) + Ai(t)76(t) + \mil{t) - a^it) - d,{t)- 
-B{t)i'i{t) - id2{t)lb{t) + Xl{t)ib{t) - id2{t)l2{t)l3{t), 
(Af (t) + Xlit)) 73 ( 075(0 - 5(075(0 + ^ 2 ( 072 ( 0 - 

-2z/ 32(074(073(0 + (cri( 0 -^i (0 + o’2(0^2(t))73(0*- 
-iA{t) - i/?2(072(075(0. 

-/3i(073(0 - *^2(075(073(0 - 2B(t)76(0- 
- 2 *^ 2 ( 072 ( 075(0 + 2(A?(0 + Ai(0)73(076(0. 
i(A?(0 + Ai(0)7l(0 - 2 *^ 2 ( 074 ( 075(0 + 202 ( 074 ( 0 - 

- |cri(t) + (cri(0Ai(t) + CT2(0A2(t))75(00 
-Pi{t)lb{t) - id2{t)iUt) + ^2(075(0 - *ai( 0 - 
-4*/32 (074(^)75(0 + 2(Af(0 + Ai(t))75(076(0 + 
+2*(cri(0Ai(0 + 0'2(0^2(0)75(0. 

-2/3i(076(0 - 2*/32(t)75(076(0 + 2(A?(0 + Ai(0)76(0. 


with initial data 

7i(0) = 72(0) = -*/o, 73 ( 0 ) = ^, 

74(0) = 0 , 75(0) = 0 , 76(0) = -^. 


(4.13) 


(4.14) 


(4.14) explicitly then we find P{t,^,x) 
Further we apply the inverse Fourier 


If we succeed to solve the problem (4.13) 
substituting ij{t),j = 1,...,6, into ( |4.12 1. 
transform to find the function P(t,f,x) and from (4.4) we obtain after integration 
f{t,x) under given initial data (4.11). 

Let us we assume that in ( |4.8[ ) 

9{x) = ^Xl-L,L]ix) = ^, xe[-L,L]. 


This choice corresponds to the initial uniform distribution of the stochastic variable 
X on the segment [—L,L]. Here f{t,x) should be read as: 


f{t,x) 


lim 

L —^ Aoo 


Ii-l.l] fPitJ,x)df 
/[-L.L] ^(^7, a^)*?/ ■ 


(4.15) 
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Consequently, 


v(t, x) 


lim 
L —^ “l“Oo 


f[-L,L] PP{tJ,x)df 


P{t,x). 


(4.16) 


Therefore the problem of finding P{t, x) is reduced to the solution of ODE system 
( |4.13[ ) with initial data 7 i( 0 ) = 0, 72 ( 0 ) = -ifo, 73 ( 0 ) = 0 , 74 ( 0 ) = 0, 75 ( 0 ) = 
0 : 75 ( 0 ) = 0 . 

Below we consider a special cases of system (4.6) arising from economic applica¬ 
tions where function f(t,x) can be obtained in an explicit form. 


4.2. Approach 2: a representation in terms of the Fonrier transform. We 

denote by P{t, /i, ^) the Fourier transform in variables (/, x) of the function P{t, /, x) 
being the solution of the problem (4.2), (4.3). Let us assume that P(t,0,^) and 


d^P{t,0,^) are functions decreasing at infinity with respect to ^ faster then any 


power of it. Then f{t,x), defined in (4.4) can be obtained as follows: 




iF^pd^P[t,Q,0]{t,x) 
p-Mp(t,o,e)](t,x) ^ 


t > 0, X £ 


(4.17) 


and ^ 


Hereinafter we denote by F~^ 
and correspondingly, let (•,-)m 
the variable /i. Here (e*^-^, 1)/ means lim (e 

L—>oo 

indicator of the set D and uJe{f) is the standard mollifier. The proof of (4.17) is an 
exercise in the harmonic analysis [30]. 

Namely, the denominator of \AA\ is 


the inverse Fourier transforms in variables 
be the action distribution on a trial function of 
^eU) * X[-L.L])/, where xa is the 


P{tJ,x)df= / Fp[Fp[P{t,^l,0]]df = 


= y^F-i[P(t,0,e)]. 

Analogously we compute the numerator: 

J fP{t,f,x)df = J fF-^[F-^[P{t,^,,0]]df = 

R R 


= f: 


,d^P{t,fi,0) ]=iV^F^^[d^P{t,0,0]- 


= V^iF^^ 


The conditional variance of P at a given value of X defined by formula (4.5) can 
be represented in terms of the Fourier transform of the joint distribution density 
P{t,f,x) as follows: 


j{t,x) = 


(P^-^ [d^P{t,0,i)]? - [c)gP(t,0,O]P^-^ [P{t,0,O] 

{F^-pp{tAmr 


{t,x). (4.18) 


We will apply the formulae to the case when the X factor volatility is proportional 
to the square root of the factor. Such model falls into the category of affine models 
[T4| therefore the Fokker-Planck equation is integrated in quadratures. 
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Remark 4.1. Affine models are popular in financial mathematics since many prob¬ 
lems can be solved analytically in their frame. In particular, affine models include 
Merton, Vasicek, Cox-Ingersoll-Ross interest rate models (see m. [H], ED). 


5. A PROBLEM OF THE PORTFOLIO SELECTION AT A FIXED TIME 

5.1. The problem statement. Let us recall that the optimal strategy by Bielecki 
and Pliska corresponds to the infinite time horizon, we are going to demonstrate 
another strategy which can be used by investor to manage portfolio and maximize 
return at any fixed time moment. 

We consider the market model of the security prices and factors (2.1), ( |2.2[ ) and 
investment process (2.3) defined in the Bielecki and Pliska model. Denoting F{t) = 
InP(f) and using Ito formula |32j we derive following equation from (2.3): 


dF{t) = 


m+n 


ito ^ lit 11 TTh 71/ 

'y ] {hiAi — —Aj y ] cr^i^) + y ( hi y ( aipXpit) 


m m+n 




i—1 p—1 




E atkdWkit), 

2 1 A) 

F( 0 ) =lnP( 0 ) = /. 


(5.1) 


We define a functional Qj{t, x; h) analogous to the first two elements of the Taylor 
series of Qe{t) about 0 = 0 in the Bielecki and Pliska model, that is 

Qj{t,x]h) = f{t,x;h) - jv{t,x;h), x = {xi, ...,Xn), (5.2) 


where 7 = | > 0 is a risk sensitive parameter analogous to 9 in the Bielecki and 
Pliska model, f(t,x]h) and v{t,x;h) are conditional expectation and conditional 
variance of stochastic variable F{t) with given values of Xi{t) = a;i, ...,X„(t) = a;„. 
Then we solve the following problem: 


to find max Q^(t,x;h), x = (xi, ...,x„), over the class of ad- 

missible investment strategies h (see Definition 5.1), with given val¬ 
ues of the factors Xi{t) = xi ,..., A„(t) = x„ at a given moment of 
time t. 


Definition 5.1. A strategy + is called optimal strategy if it gives a maximum of 
the functional Q-yft, x; h) with given values of the factors Xi{t) = xi,..., A„(t) = x„ 
at a given time t. 


Once we find the maximum over the class of denoted strategies then we find 
the strategy which provides the maximum portfolio return with regard to loss of 
random nature described by the variance. Changing the value of parameter 7 , we 
can overstate or understate the role of randomness, or do not take into account the 
randomness at all, setting 7 = 0. The model can be interpreted in the following 
way. Let us assume the investor be going to allocate the initial capital between 
assets Si, i = l,...,m, with the prices depending on a set of exogenous economic 
factor Xj, j = l,...,n. The prices of assets and values of factors obey equations 
(2.1), (2.2). The investor solves a dynamic asset management problem featuring a 
risk sensitive optimality criterion. Let us assume that the investor knows an explicit 
values of factors in a fixed moment of time. Thus, the investor has to find the 
optimal portfolio taking into account a new information on the factors, such that 
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the model is flexible and can be actualized within all the time of investment. The 
model refers to the tactical asset allocation (e.g. [55]1. 


5.2. The algorithm of solution. Let us give a outline of solution to the optimiza¬ 
tion problem for a model with one factor, that is we consider system (5.1), (2.2) for 
n = 1 (here Xi{t) = X(t)): 


dF{t) = 


cr^fc) ^ hia,X{t) 


m +1 


k^l 


djt-\- 


m m+1 


(5.3) 


a^kdWu{t). = 


2 = 1 


m+1 


dX{t) = {B + l3X{t))dt -f XkdWkit), X{0) = x. 


(5.4) 


k=l 


We can use formulae (4.4), (4.5), (4.15), (4.16) or (4.17), (4.18) to find f{t,x;h) 
and u(t, x; h). 

Then we can write Q^it, x; h) as a quadratic function with respect toh — (hi,..., 
Below we will write this function explicitly for several important cases. 

m 

To find a conditional extremum of Q-y(t, x, h) with the constraint — 1 = 0 


the Lagrange method can be applied. The Lagrange function is 


i=l 


L{h,i) = Q^{t,x\h) + 

2=1 

m m 

= Y Kij{t,x)h^hj + y^(iLi(t,a;) -f Ko{t,x) - 

2 ,^ = 1 2=1 

where Kij, Ki, Kq are functions of t,x and coefficients +, B, /3, Afc, i = 
1,..., m, k = 1, ..., m -I- 1. 

We get a system of m -I- 1 equations by equating to zero partial derivatives with 
respect to hi,^ of the Lagrange function L{h,^): 


dL{h,0 

dhi 

dL{h,0 


Y, i.Kij{t,x) + Kji{t,x)) hj + Ki{t,x) - 1 -^ = 0, 
i=i 

m 

^ h, -1 = 0 . 

i=l 


This is a nonhomogeneous system of linear algebraic equations with respect to vari¬ 
ables hi,..., hm, C The unknown hi,..., h^, ^ can be found uniquely provided the 
determinant of the system does not vanish. If lim Q^it, x; h) = —oo and Q^{t, x; h) 

|ii|—>00 

is continuous function in h, then the point hi,..., hm is a unique maximum. 


Remark 5.1. Our considerations hold for 7 > 


1 

2 ' 


Remark 5.2. We restrict ourselves by consideration of the model with one factor 
X{t), since our main goal is a study of influence of such factor as the interest rate. 
Nevertheless, the results can be extended to the case of vectorial equation with n 
components for the factor process X{t). Here f{t,x) and v{t,x) are functions of 
time and n spatial variables. The factor process can have correlated components. 
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Further we find an explicit optimal strategy of investment for the case of portfolio 
consisting of two assets, depending on one market factor, the bank interest rate. For 
the interest rate we choose first the Vasicek model and then the Cox-Ingersoll-Ross 
model. 


6 . The linear interest rate (the Vasicek model) 


6.1. An example of portfolio consisting of two assets. Let us write (5.3), (5.4) 
in a general way: 

dF = (A + aX)dt + {a,dW), 

dX = {B + (3X)dt + {X,dW), (6.1) 

F{0) = /, X{0)=x, 

where 

m ^ m+1 m 

A = ^ ^ V ^ik)^ ^ ~ ^ ^ A = (Ai, Atti-I-i), 




— ^ hjCTil-j ^ ^ hjCTj^ 


m+1 J 5 




2=1 


( 6 . 2 ) 

W = (IFi(t),..., Wm+i{t)) is a (to + 1)- dimensional Brownian motion. Recall that 
/3 < 0. 


Equation (6.1) is a particular case of (4.6). Therefore we can use the result of 

Secim 


(6.3) 


Equation ( |4.7| ) is 
dP{t,f,x) 


dt 




df 

1 d^P{t,f,x) 1 d^P(t,f,x) 




dp 


-E 2 - 

2 


dx'^ 


3x 

PP{tJ,x) 

* dfdx ' 


where 


m m+1 

m+1 

E2 = AA^ = ^ Xl 


2^1 k—1 

m+1 m 

fc=l 

(6.4) 

crA'^ = ^ Xky^pjCTik- 




k—1 i—1 


The initial conditions are 


PiOJ,x) = S{f-fo)gix). 


To solve this equation we use the first approach from Sec|^ nevertheless, the 
second approach can be applied, too. We will show how this second approach works 
in Secj^on the example of the Cox-Ingersoll-Ross interest rate. 

6.1.1. Gaussian initial distribution of the factor. To get explicit formulae we con- 

e ^ 

sider Gaussian initial distribution of the random value X, that is g{x) = - -j =—, 

v27rs 

where Xq G K is the mean value of X at initial moment of time, the constant s^, 
s G K+ is the variance. The limit case s —?► 0 corresponds to the factor that equals 
xo initially. Thus, 

_ 

P{0,f,x) = 5{f-fo f " . 

y/ZTTS 


(6.5) 
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The Fourier transform with respect to / maps (6.3), (6.5) into 


dP{t, x) 


+ ax)P{t, fi, x) — {B + fix )-— ^ — PP{t, x) — 


1 2 " / \ 1 d'^Pit, /i, x) 

-SiM P{t,fi,x) + -Y, 2 - - ^ - 


dx 

+ i/iSa 


dP{t, x) 

dx ^ 


( 6 . 6 ) 


P{0,fj,,x) = 


_ (x-xq)^ 
g 2 s^ 


■\/^S 


(6.7) 


The anzats for the solution to (6.6),(6.7| is 


P{t,fi,x) = 


o7l(*)+72(t)At+73(t)2:+74(i)l‘^+75(i)2'M+76(t)a:^ 


•\/^S 


( 6 . 8 ) 


We substitute (6.8) into (6.6), (6.7), then equating the coefficients at the same 
powers of fj, and a;, we get a system for 7 j(t), j = 1,..., 6, a particular case of (4.13) 


d'^2it) 

dt 

dizjt) 

dlkt) 

dt 

djait) 

dt 


'P‘2l&{t) + 5^271(1) - /3 - ^73(1), 
S273(075(i) - B^^{t) -iA + iT.z^z{t), 
2^273(i)76(i) - 
— ^Ei + ^£375(1) + 4 E 27 |(t), 

2*^376(t) -ia + 2 ^ 275 ( 1 ) 76 (t) - dibit), 
-2/376(1) + 2 E 27 |( 1 ), 


(6.9) 


subject to initial data 


7i(0) = -^,72(0) = -*/o,73(0) = f^, 

74(0) = 0,75(0) =0,76(0) = -2^. 


( 6 . 10 ) 


Solving problem (6.9), (6.10), we get 7 j(l),j = 1,...,6, explicitly: 
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71W = - 


72 (t) = 


73 (i) 

74 (i) 


75 (i) 

76 (i) 


+ ^E2 In ( _E2+2e24t^/3s2+e2^<ts2 ) 

2 ^(-E2 + 2e2/3*/3s2 + e2/3‘S2) 

(S2 + 2/3 s 2) In ( -S2+2e^^yj^+e^'^*S2 ) g2/3i 

2 (-E 2 + 2e2/3‘/3s2 + e2/5*S2) 

2B{B + xo/3)e'5* + (B + 

^ /3(-E2 + 2e2/3i^s2 + e2/34E2) ’ 

, (2Ba + a:o/3a + PBat - l3^At)T.2 - 25^3/3 
-^Jo - ^ - 


—i 


-i- 


( 2 /?s 2 + S2)/32e2/3‘ - E2/32 

e/5«((-4Ba - 2xoPa)^2 + + 2 /? 2 a;o)E 3 - 2s'^j3Ba) 

( 2 ^s 2 + E2)/32e2/3* - E 2/?2 

e^^*((i?a(2 - /3t) + /322lt + aio/3a)E2 - 2f3{B + /3a;o)E3) 




( 2 ^s 2 + E2)/32e2/3‘ - E 2/32 
e^^*-{—2Bj3‘^s‘^at + 2s'^j3Ba + 2s'^/3^At) 

'^'^{2Ps^"+t2)We^^*-~^2W ’ 

2(-B + e^‘a;o/3 + 6/^*5) 

-S 2 + 2e2^**/^s2+ e2^*E2’ 

Eli —2(E2a — E3/3)2 + Yi2aj3{2'Ts^j5t — as^ — Yi2at) 

~Y~ ^ (2/3s2 + E2)/33e2/3i - E2^3 

4e^‘(E2a — E3/3)(E2a — E 3/3 + a/3s^) 

^ (2/3s2 + E2)/33e2/3* - E 2/33 

g2/3t ((2 s 2^2^ ^ _ 3^2)^ _ 2E2)E2a2) 

(2/3s2 + E2)^3e2/3‘ - E 2/33 

^ e2/5*2 {{-2s^pH + (2s2 - E2t)/3 + 2E2)E3a - 2E2/?) 

( 2 ^s 2 + E2)/32e2/5‘ - E 2/32 ’ 

2 i{aj 5 s^ + q;E2 — E 3 / 3 )e^‘ — za(E2 + 2 ^s^)e^^* — 2E3/3 + q;E2 


+ 


/? 


((2/3s2 + E2)e2/5* - E2)/3 


-E 2 + 2e2/3‘/3s2 + e2/5*E2 


( 6 . 11 ) 

We substitute the expressions for = l,-.., 6 , in ( 6 . 8 ) and find P{t,fi,x) 

explicitly. This expression is cumbersome, we do not write this. The inverse Fourier 
transform gives 


P{tJ,x) = ^= P{t,fi,x)P-f>^dfi = 

V J —00 


P'2gCi{tJ,x) 

yPrsC2{t) 

with 

C2(t) = (|2/34l]its2 + (SEgats^ + EitE2)/33+ 

+(—4s^E2Q;^t + 4E| — 8 s^q;E 3 + 4aE3tE2)/3^+ 

+(-8aE2E3 - 2a2E2t + 6a^s^E2)P + 

+ ((-8E2 + 8s2q-S 3)/32 + (16aE2E3 - ia^s^Y.2)P - - 

-EitE2^3 ^ (_4aE3tE2 + 4E2)/32+ 

+ (2a2E2t — 8 q:E 2E3 + 2a^s^E2)/3 + 4 q;^E2, 


Ci{t,f,x) is a function of variables t,f,x and parameters /o, Xq, s, A, a, B, (3, Ei, 
E 2 , E 3 . The multiple integration by parts was performed by means of the computer 
algebra system MAPLE. 
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The integral converges under condition 


:= ( - 4a^El - 4 S §^2 + 8aS2S3/?+ 


U := 


+ (E 2 + 2/3s2)(2a2E2 - da/SEg - 


= 2 / 3 t _ 


(E 3/3 - E 2 a)^ - s'^al3{T,3/3 - T,2a) 




+8 

+ (Ei/32 - 2a2E2 + 4aE3^)E2/3t - 4Ef/32 - 4(E3/3 - E2a)2- 

-1 


-2a^s‘^^2P) ((2/3s2 + E2)e2/3* - E 2 ) > 0 . 


( 6 . 12 ) 


Since C/(|t=o = —Ei/3^ > 0, then there exists t* > 0, such that for all t G (0,t*) 
this condition holds. Nevertheless, lim U = 00 , the sign of the infinity is equal to 

t—¥00 


the sign of (—Ei/3^ 
satisfied for large t. 
Let us denote 


4aE3/3^ + 2E2a^/3), therefore the condition (6.12) can be not 


Us^ := lim U = 

= 4aE3/3 - 3E2a2 + {2l3T.2a^ - Ei^^ - 4/3^T,3a)t+ 

+(-4aE3/3 + 4E2a2)e-'5* - a'^T,2e-^^\ 

UsQ ■■= lim U = 

-4E|/32 + 8E2a/3E3 - 4a^Sl (EiEa/S^ + 4aS3^2/3^ - 2a'^El/3)t 
~ (e2/3i _ i)E 2 (e2/3t _ 1)E2 

4(E2^2 _ 2E2a/3E3 + a^Y;l)e^P* i2a^^2P - Ei/?^ - 4aT.3P‘^)e‘^^H 
(e2/3t_i)E2 

(8Ei/32 - 16E2a/3E3 + 8a2E2)e/5‘ 

^ (-l + e2/5‘)E2 ■ 


Therefore for large s condition (6.12) is not satisfied as t —?► 00 , since lim = — 00 . 

^ -' t —>-00 °° 

For small s the condition (6.12) is not satisfied for all t > 0 if —Ei/3^ — 4/3^E3a + 
2a2S2/3 < 0 . 


Then we substitute P{t, f, x) in (4.41,(4.5), integrate and get the expectation and 
dispersion: 


f{t,x) = —(^ {2j3'^xa + 2Bali)e^^+ 

+ {-2l3^xa + (-27l< - 2 /o)/33 + 213'^Bat - 2Bal3)e^l^^ 


+2 


(E 2 Q: — E3/3)/3a; — /3^XoE3 + (E2a:oQ; — 2BYj3)f3 + 2'Bi2Ba 


=/ 3 i_ 


- PxaT,2 + ( 2 X 0 E 3 - 24tE2 - /oE2)/3^+ 

+ ((BtE 2 - E 2 a:o)a + 2 BE 3 )^ - 2 E 2 Ba] 

+ (—E 2 / 3 Q; + 2Ti3P^)x + (AtE2 + /oE2),0^+ 

+ ((—E2a;o — Btl]2)a + 2 BY, 3 )f 3 — 2Yi2Ba^ x 

X (/ 32 (-E 2 + 2 e 2 / 5 */ 3 s 2 + e 2 / 5 ‘E 2 )) 


(6.13) 
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~j{t, x) = (^ 8(/3^S3a — E2a^/3)e^* + 2Il2a'^/3+ 


+ {2Eitl3^ + - 4(2E3a + + QE2a‘^ 3^ + 


+8 


2E2aE3P - Elo? - El0^ 


=/3t , 


4(S3affi2 + E§)/32- 


-2(4aE3S2 + ElaH)l3 + + Ei0^E2 

-4(E3atE2 - E§)/32 - 2(4aE3E2 - ElaH)/3+ 


=2/3t_ 


(6.14) 


+4Eia2 _ Eit/33E2) 


-E, 


■ 2e2^*/3s2 ■ 


■e2/5*E2) 


-1 


From (6.13), (6.14), (6.2), (6.4) we get an explicit formula for Q^{t,x]h). 

Let us consider the optimal strategy on the example of two assets, depending on 
one market factor, the linear interest rate. It this case it is convenient to write the 
strategy in the following way:(/ii,/ 12 ) = (h, 1 — h). Then from (6.13), (6.14), (6.2) 
and (6.4), where m = 2, hi = h, h 2 = 1 — h, we get: 


Qy{t, x] h) = K2h^ + Kih + Kq, 


(6.15) 


where 


K 2 = 


(8[(/3E2s 2 + E2)(ai - a2f - (^d^s^ + 2pE2)S2(,ai - ^2) + 51/3^ 

((—4 + 2l3t)E2 + (—6/3 + 4/32t)s2E2)(ai — 0 : 2 )^ — P^t(E2 + 2/3s‘^)S3+ 


+4 ((2 - /3f)/3E2 + 2(1 - /3t)/32s2) ^^(ai - a2) - 


ye 


2/3t_ 


I3H{E2 + 2I3s^)Si 


=2/3t , 


- 2((2 + pt)El + /3E2s2)(ai - a2)^+ 
+4{pH + 2 / 3 ) 52 E 2 (ai - ^ 2 ) - 4^|^2 ^ ^ _ pHSiE 2 'j x 

x(/33(2e2'3*/3s2 + (-l + e2/3*)E2)) \ 


(6.16) 


Ki = (^- 8 / 32(1 - eP^f-iS3S2 + [ - 8/3a2(2E2 + /3s2)e/3* + 

(4/3a2(2 — f3t)E2 + 8 / 320 : 2(1 — I3t)s'^) + 4,002(2 + /3t)E2 

2{{xo + x)p + 2B)^2g/3t _ 2{^pxo + B)f3^e^^^ - 2l3‘^{(3x + B) 


1 S 2 + 
52 + 


{2I5^sH + /33tE2)e2/3‘ - j3HE2 


Sa 


- 8(;3s 2 + 2E2)/3e^‘+ 


+(4/3(-/3< + 2)E2 + 8 / 32(1 - /3t)s2)e2/3* + 4(/3t + 2)/3E2 (oi - 02 ) 7 ^ 5 + 


{- 2 PHE 2 - + 2 PHE 2 


'jSq + 4a2E2 


4(E2 + /3s2)e/5*+ 


+((/3f - 2)E2 + /3(2^t - 3)s2)e2/3‘ - E2(2 + /3t) - (3s^ (oi - 02 ) 7 - 
-2 ((/3(x + xo) + 2B)PE2 + iPx + B)P^s^) (oi - 02 ) 6 ^ 5 *+ 

+ [{2B - Bpt + (a: + xo)j3)l3E2 + 2{B - B(3t + Px)P‘^s‘^) (oi - a2)e2^*+ 
+(i?/3t + (a; + xq)P + 2 i 3 )/ 3 E 2 (oi — 02 ) + {Ai — A 2 ){E 2 {&^^^ — 1) + 2/3s2)^ x 
x(/33((E2(e2'3‘-l) + 2/3s2))“\ 

(6.17) 
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Kq — ^ 2 Yl 2 Bo £2 + (— 2 S^B + a2^2(xo + 2; + Bt))j 3 -{- 

+{-285X0 + S 2 (/o - 1^4 + ^) + 2s^a2ix - Bt))P^e^P^- 
—2{a28l2ixo + a;) + Ba2s'^ — 2BS5)j3e^*+ 

+2{S5{xo + x) - s2xa2)/3^e^‘ + 2s^(/o + ^ - 5S'4)e^^‘/3^+ (6.18) 

+ (S 2 (i ELi ^ - /o) - 2 x55 )/3^ + 

+ (a 2 S 2 (—.Bt + xo + x) + 2 Ba 2 S^ — 2i?S'5)e^^*/3+ 

+2Ba28,2{e^l^* - ({2e^l^^Ps^ + ^2(6^^* - l))/3^) \ 

=-IELi'^ 1*+ 0-2*: Ba =Ei=i('^i*-<^2*)A*, 53 = Ei=i(^n 

^3 2 a v^3 _ \ o v^3 


*5'4 — Ei=l ' 


5*5 — 1 Sq — 1 *^22(^11 ^^22)’ 


Since ^2 = 0 at 2 = 0 and ^\t=o = -7ELi(^i* “ '^ 2 i)^ - 5 ELi('^ii + ^^iz) < 0, 
then there exists t* > 0 such that for Vt S (0,t*) a unique point of extremum 
Hj = is a maximum point of Qj{t,x; h) with respect to h. 

Thus, knowing the value of the interest rate at any moment of time t S (0, t*) we 
can maximize the income investing a part in the first asset and the rest (1 — H^) 
in the second one. 

In the limit cases s —^ 0 and s —^ 00 the expression Q^{t, x; h) becomes simpler: 

(3 + e-2/3* - 2j3t - 4e-^‘)S27a2 


lim Qj{t,x;h) = fo + tA —— 


13^ 


4(/3t — 1 + e ((/?x + B)(e — 1) +/3Bt) 


\iu\Q^{t,x-,h) = At + fo—^tTii + 


j3'^ /32 

{2{pt-2)e^^+pt +2)^210^ 


s—>-0 


-4 




((/3t-2)eP‘ + /3t + 2)S37a 
/32(e/5* + 1 ) ^ 

{{{—Bt + Xo + x)f3 + 2B)e^* — (xq + Bt + x)j3 — 2B)a 
/3'^{e^* _l_ 

4(1 - e^‘) 7 E| 2i-{B + /lxo)e^* + /3x + 5)113 


/3(e/5‘ + 1 )E 2 


f3{eP* + 1)E2 


where A,a and Ei,S 2 ,S 3 are given by ( 6 . 2 ) and (6.4), respectively. 

6.1.2. Uniform initial distribution of factor. We consider a random value X dis¬ 
tributed uniformly on the segment [—L,L], that is g{x) = i^X[-l,l]{x). Then 

n/n J- ^ - fo)Xl-L,L]{x) S{f - fo) r rl i n^ 

5(0, /, x) = 2L ^ -TL-’ 


After the Fourier transform with respect to / the equation (6.19) takes the form 
P(0,M,x) = xG[-L,L]. (6.20) 


For the solution of problem (6.3), (6.20) we use the anzats: 

g7l(i)+72(t)M+73(t)2:+74(t)A‘^+75(i)a:'M+76(t)a:^ 


P{t,n,x) = 


2LV^ 


( 6 . 21 ) 


We have to solve (6.9) with initial conditions 71 (0) = 0 , 72 ( 0 ) = — 1 / 0 , 73 ( 0 ) = 
0 , 74 ( 0 ) = 0 , 75 ( 0 ) = 0 , 70 ( 0 ) = 0. Thus, 7 j(t),j = 1, •.•, 6 , can be found explicitly: 

, , „ , , aBit aBie~^* aBi ^ „ 

li{t) = -!3t, 72 ( 2 ) = ^—-k— - -^-Ati-foi, 73(t) = 0, 
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_ ( —4Xl2a^+4aS3/5) . a^S2ea:p( —2/3i) . (4aiS3/3 —2Sit/5^+3E2a^—2Q^tS2/3—4aS3/3) 


/,x _ ^-42.2a +4:a^3P) I 
74l(-j — 4/33ea;p(-/3t) 


4/33 

ai 


+ 


4/33 


75(0 = -j{e - 1), 76(i) = 0. 


We substitute these expressions to (6.21) and find P{t,^,x). The inverse Fourier 
transform gives us P{t,f,x). 

The following restrictions have to be imposed to guarantee the convergence of 
integrals for t > 0 : 


- 4(a2S2 - a/3E3)e-^‘- 

- a^T.2{2l3t - 3) - 4a^S3(l - pt)j > 0, 

- + 2a^l3Y,2)t + 4aS3/3- 

-3a2S2)e2/5* - (-4a^E3 + 4a2S2)e^‘ + a‘^T .2 > 0. 


( 6 . 22 ) 


(6.23) 


One can show that inequality (6.22) takes place for any parameters, whereas for 
(6.23) is true only for j3 < 0. 

Thus from (4.15), (4.16) we get the following values for the conditional expectation 
and variance: 


j., . iPx + B)ae {Ba - Aj3)t {fix + B)a 

f{t,x) = -- - -+--+/o, 


(6.24) 


v{t,x) = — 


(-4S2a2 + 4aE3^)e-^‘ E2a2e-2/3‘ 

2/3'3 ■■ 2)P 

{4:a/3Hj:3 - 2Eit^3 35^20,2 _ 2a^ptY.2 - 4aE3/3) 


where A, a, Ei, E 2 , E 3 are given in (6.2), (6.4). 

For the case of two assets depending on one factor we get 

Q.yit, x; h) = K 2 h'^ + Kih + Kq, 


(6.25) 


where h and {1 — h) are the proportions of capital invested in the first and second 
assets, respectively , where K 2 , Ki,Ko are functions, which expressions are cumber¬ 
some, nevertheless, the dependence on time is only exponential or polynomial. Since 
K 2 = —^ {<^ ii + (^ 2 i) ~ iv{t, X] h) < 0 , then the optimal strategy in the sense 

of Definition 5.1 is In the next section we find the explicit optimal 

strategy for the classical example of portfolio containing of two assets one of which 
is a bank account. 


Remark 6.1. We performed a series of numerical experiments and showed that for 
a real market parameters the difference in optimal strategies depends very weakly 
on initial distribution of the factor. Namely, for the Gaussian distribution with 
s > 0.001 the result is very similar to the limiting case of uniform initial distribution. 

6.2. Comparing with the Bielecki and Pliska strategy. T.Bielecki and S.Pliska 
in their works considered a classical example of portfolio consisting of two assets, 
where one asset is a bank account and a factor is the interest rate. The formula for 
the Bielecki and Pliska optimal strategy Hg and the maximal value of the functional 
p{6) is written out in Secj^ To compare our strategy with the Bielecki-Pliska one 
we also consider this classical example. 

Thus, let the assets of the portfolio obey SDEs: 

dSi{t) 


{Ai + aiR{t))dt + aidW\{t), iS'i(O) = s > 0, 
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dS2it) 


= R{t)dt, 52(0) = 1, 


S2{t) 

the dynamics of the interest rate R is 

dR{t) = {B + l3R{t))dt + \dW 2 {t), i?(0) = r 

(the Vasicek model). Here Hi, oi, /3, cti, A are given constants, moreover B > 
0,/3 < 0, and Wi,W 2 are independent Brownian motions. 

The equation for the capital of investor is the following: 

dV(t) 


V{t) 


= [hiAi + {hiai + h 2 )R{t)] dt + hiaidWi{t), V{t) = v > Q. 


Since we consider the portfolio consisting of two asset, the we denote hi = h the 
proportion of capital invested in the risky asset and ^2 = 1 — h the share invested 
in the bank account. 

If lnl/(t) = F{t), then 


dF{t) = 


hAi -+ {hai + 1 — h)R{t) 


dt + haidWi{t). 


We consider the case of uniform initial distribution of the interest rate R. Here 
(6.24) (6.25) have the form 

f{t,r) = M 2 h?{t) 


where 

M 2 

Mo 

Li 

Then 


Mih(t) + Mo, v{t,r) = L 2 h?{t) + Lih{t) + Lq, (6.26) 

alt (oi - l)(/3r+ B)(1 - e"^‘) (H(ai - 1) -/3Hi)t 

- j- -:- -p -■ 

(/3r + H)(l - e ^‘) Bt 

-^(ai - l)<(>(i), Lq = Ht) = (e"^^‘ - 46"^^* - 2/3t + 3). 


L2 = lf(t){t) + alt, 


Q-y{t, r; h) = {M 2 - "fL 2 )h'^{t) + (Mi - 'yLi)h{t) + Mq - yLp. 


(6.27) 


o-Pt _ 


1)" 2 

- \- ai > 0 , then the coefficient of 


c- T rn\ n dL2it) ^^{^1 - l)^(e 

Since L 2 ( 0 ) = 0, —^ 

the leading term of the quadratic with respect to h function Qj{t, r) is negative and 
the unique point of maximum is 

^ Ml -^Li 

^ alt + 27L2 

Thus, at any moment of time the investor, knowing a current interest rate, can 
maximize his income investing a proportion of capital to the risky asset and the 
rest 1 — to the bank account. 


It follows from (6.26) that as t —>■ 00 the conditional expectation f{t,r) and the 
conditional variance v{t,r) increase as and respectively (we recall that 

/3 < 0). Introducing the risk coefficient, we describe the subjective influence of 
randomness to the expected mean income of portfolio. 

Fig§ shows the dependence f{t,r-,Hj) on v{t,r-Hj) (the effective frontier) at 
different moments of time for different values of the parameter of risk 7 and given 
other parameters of model Hi = 0.15, oi = —l,ai — 0.2, B = 0.05,/3 = —1,A = 
0 . 02 , r = 0 . 01 , /o = 1 (the values of parameters are chosen as in example from [^). 

Let us compare the conditional expectation of the portfolio at a fixed value of 
factor for two strategies under consideration. We substitute (6.28) and (3.7) in the 
formula (6.26) for f{t,r) and after computations we get the following proposition. 




















22 


G. S. KAMBARBAEVA AND O. S. ROZANOVA 



Figure 1 . as a function of v{t,r;H^) at fixed moment 

of time and different 7. 


Proposition 6.1. For 7 = 0 , d = 0, t>0 the following inequality holds: 


{f{t,r;H^) - f{t,r-Hg))\^=o^e=o = 


2a1 


> 0 . 


Proposition 6.2. There exists t* > 0 such that for all t S (0,t*) the following state¬ 
ment hold: 

( 1 ) if ai ^ 1 and the factor satisfies the condition 

then there exists 7 * > 0 such that for all 7 G ( 0 , 7 ,) 
f{t,r;H^) - f{t,r;Hg) > 0 ; 

(2) if ai = 1 and Ai > 0, then f{t,r;Hj) — f{t,r]He) > 0 for all 7 > 0. 

Proof. We consider the difference f{t,r;Hj) — f{t,r]Hg) as a function of t at 
other parameters fixed. We denote this difference as q{t). The computations show 
that ( 7 '( 0 ) = 0 and 

g"(0) = (8r-VA(ai - 1)3 - 2((1 + 27)7/3air2 + 

-|-((1 -I- 2j)'yB(7i — 87 ^AAi)r)(Q;i — l)^ — 

—2^(— 2(27 -I- l) 7 ^AtTi + (1 + 27)7Ai/3tTi)r-|- 
-|-(1 -I- 2^)^AiBai — 47^AAJ^ (oi — 1)-|- 
-h 4(27 -h l)j'^XajAi){crl{l -h 67 -f 127 ^ -p 873 ))-!. 

The denominator of this expression is positive. One can see that if oi = 1 and 
Ai > 0, then (;"(0) > 0. This proves the second part of the proposition. 

If oi ^ 1 , than we expand the numerator (;"( 0 ) = g"( 0 , 7 ) in series with respect 
to 7 about zero: 

g"(0,7) = -2/3cri(ai - l f{r+ ^){r + ^^ )7 + Q(7^)- 


This expression is positive if the condition (6.29) holds. □ 
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Figj^ shows the graphs of conditional expectations for two strategies for the same 
parameters of model as in Figj^ Here 4j = 0 = 0.1. This example shows that 
our strategy for realistic values of parameter gives greater expectation of portfolio’s 
long-run expected growth rate till the moment t* and this situation can hold within 
several years. 



Figure 2. The conditional expectation of portfolio’s long-run ex¬ 
pected growth rate h) for our strategy (solid line) and in 

the case of the strategy of Bielecki and Pliska Hq (dashed line). 


Taking into account the principle of constructing the function Q^{t^ r; h) it makes 
sense to compare two strategies for 7 = 9/4 and small 6. First we compare the 

results as t —>■ 00 . At any fixed x we get lim H^{t) =- - in the case ai 7 ^ 1 


ai — 1 


and lim H^(t) = in the case ai = 1. Thus, the limit Hj(t) is discontinuous as 

t —^CXD OX 2 

a function of ai. We introduce the following denotation: 


P{l) ■= lim 

t—>oo 


t 


Computations show that 


p{l) = t=FT-(7+2)(FI^> 

P(7) = 5 


B 

/3’ 


tti = 1. 


The value ^( 7 ) corresponds to the expected rate of growth of capital at infinity, it 
is analogous to p{6) in the model of Bielecki and Pliska, see( |3.8[ ). Figi shows this 
function at the same values of parameter as at FigQ 9 varies from 0 to 1, 7 
After the limit pass as 0 —>■ 0 we get from (|3.9[) 


0 

4- 


B 


1 


B 


l.m „ = +^(.4,--(01-1)7- 


A"(ai - 1)= 
4cr^/3 


(6.30) 


If we set directly 7 = 0 in ^ then perform the limit pass as t —t 00 , we 


get (6.30) without the last term, containing A. 
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Figure 3 . The expected rate of growth of capital at infinity: ^(y) 
for our strategy, p{9) for the Bielecki and Pliska strategy, 9 = dy. 


6.3. Asymptotics of the proportions of capitals. In this section we study the 
asymptotics of proportions of the portfolio capital as times goes to infinity for the 
cases of two and three assets depending on one factor with uniform initial distri¬ 
bution. We recall that for the Gaussian initial distribution basically there arises 
a restriction on application of strategy for large t, therefore the computation of 
asymptotics is impossible. 


For the case of two assets depending on one factor we have the system (5.3), (5.4) 


for m = 2. The functions corresponding to the proportions of capital are found in 

SeciHa 

= (6.31) 


H2{t) = 1 + 


2^2 ’ 
Ki 
2Ko 


(6.32) 


We denote the asymptotic limits of the proportions of capitals as follows: 

HT(2)-= / 4 m : = ^lim (^^ 2 (t))■ 

The computation shows that the following proposition holds: 


Proposition 6.3. Let Hi{t) and i? 2 (t) be defined as (6.31), (6.32), then 
(1) if ai 4 0 ^ 2 , then ~ ~~ 


JLfOO _ 1 ZJOO , 

^2{2) — ^ 1 ( 2 )J 


ai — Q!2 

(2) if ai = a 2 , then H^ 2 ) = • 00 , 

_ 7Q;i((crii — (Ti2)Ai -I- ((T21 — cr22)A2 -I- ((T31 — cr32)A3) 

^ (0-4 -k cr4 + 0-13 + ^21 + <^22 + 0-13) (2/37 - 1) 


(3) if JCi = 0, then Hi{t) = = 


(2/37 ~ I)®”! + ^2 ~ ^1 


( 2 ) ( 0 - 2 + cr|) (2/37- 1 ) 

Thus, the limit values of proportions at infinity in the case of two assets depends 
only on parameters ai and 0 : 2 , if cti 4 « 2 - In the case ai = 0:2 the limit value 
depends on other parameters, too. For large t it worth to invest to an asset that is 
less depending on the factor (the respective Ui is the smallest my modulus), despite 
of trends and volatilies. 

In the case of three assets depending on one factor (m = 3, n = 1) the proportions 
of the capital can be also found according to Sec|^ Let us introduce the denotation 
of limits of proportions of the capital Hi(t),i = 1,2,3 as the time tends to infinity: 

'■= A™ ■- lim {H 2 {t)),H^^^ := lim (Hsit)). 
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The explicit formulas for asymptotic limits are the following: 

^i{3) ~ (2/?(C^3(<T2 ~ Q:i<T2) + U2{a\ — aiQ;3))7 + (~^i + ^3 ~ U^)a\+ 

+ ( — 2 I 1 + 2 I 2 —172)ci3 + (2^1—^2~^3)<T2<T3 + (^2~^3 + C^3)<T1<T2 + (~^2+^3 + C^2)<T1<T3)/ 

{2j3{{U2 + U3)a\ + {Ui-\-U3)a\ + {Ui + U2)oi\)+Aj3{—U3aia2 — U20tia3 — Uia20i3))^+ 

+{— U 2 — Uz ) oi \ + {—Ui — U ^) a \ + {—Ui — [72)0:3 + 2 ([/ 3 Q;iq ;2 + [720103 + [710203)), 

^^ 3 ) ~ (2/?([73(o3 — 0102 ) + [7i(o3 — 0203))7 + ( — 7 I 2 + 2 I 3 — [ 73)03 + 

+ (7ll— 7 I 2 —[7i)o 3 + ( —Ai+2yl2—7l3)oi03 + (7li—7l3 + [73)ai02 + ( —^1+^3 + [7 i)0203)/ 

( 2 / 3 ( (7/2 + [73)03 + ([73+772)02 + ([73 + 7/2)03) + 4/3 (—[/3O1O2 — [72O3O3 — [ 7 i 0203 )) 7 + 

+ ( —[ 72 ~C 73 )o 3 + ( — [73 — [[3)02+ (—[ 7 i — [[2) 03+2(77303 02 + 77203 03 + 7/30203)), 


-^+ 3 ) “ (2/?(C72(o 3 — O1O3) + [/l(o2 — 0203))7 + {A2 — — 7/2)03 + 

+ ( 2 li— 2 I 3 — 7/3 )o 2 + ( —^1—24.2+ 2 ^ 3)03 02 + (Al— 7 l 2 + 7 / 2 ) 0303 + ( — 211 + 212 + 773 ) 0203 )/ 
(2/3((7/2 + 773)o3 + ([7i + 772)o2 + (77i+7/2)o3)+4/3(—7/301O2 —772O1O3 —77i0203))7+ 
+(— 7/2 — 7 / 3 ) 03 + (—77i — 7 / 3 ) 02 + (—77i — 772 ) 03 + 2(7730102 + 7720103 + 7 / 10203 )), 


where 


TJ _ 2 I 2 I 2 I 2 

— ^11 + ^12 + ^13 + *^14 5 


U 2 — (721 + <^22 + ^23 + '^24 7 




31 


^2 I ^2 I ^2 

(^ 32 + “r 


^33 ^34- 


Let US note that the limit behavior depends in general case on all parameters of 
the model and this difference from the case of two assets seems strange. Nevertheless, 
if parameters Oi for a pair of assets coincide, then the situation is analogous to the 

case of two assets. For example, if 02 = 03 , then for any values of other parameters 
d 2 

= -. Moreover, and T/^g) depend on other parameters of the 

model, too. The case where all Oi are equal, is degenerate, as above: 


H 


1(3) 


= /Cl • 00 , /Cl = 


7C^i(Eti Ai - ^, 1 ) 

(2/?7-l)EkLi=i^M) 


where cr/ = i = 1,2,3, where i,j,k are all even combinations of indices 


fc=i 


(1,2,3). If/Cl = 0, then 
H,{t) = T/i^g) = 


( 2/37 - 1 )(t|(t| + ( 2 I 2 - ^i)cr| + {A 3 - Ai)a^ 


( 2/37 - 1) j=l 
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6.4. Influence of different parameters of model on the optimal strategy 
of investment for small time. It is interesting to note that for small time the 
strategy of investment depends on all parameters and differs significantly on the 
limit behavior as t —>■ oo. We show the results of computations for the case of 


two and three assets (the values of parameters are given in the tables 6.4 and 6.4 
respectively). 

As we have seen for the case of two assets (n = 2) the limit behavior of the 
proportions depends only on parameters 01,02- For small times the strategy is 
different: 


• Fig. |4] presents graphs for the proportions of capital for different (3. For 
greater (3 the function reaches its asymptotical value quicker; 

• Fig. 0 presents graphs for the proportions of capital in dependence on pa¬ 
rameter (Til. For small t an increasing of an results to a decreasing of the 
proportion of corresponding asset in the portfolio; 

• Parameter Ai influences very weakly at large times, nevertheless, for small 
t its influence is significant (see Fig. [^. 

To study the strategy of optimal investment in the case of three assets (n = 3), 
we analyze functions = 1,2,3, changing the values of parameters /3,aij,Ai, 

f = 1,2,3; j = 1,2,3,4: 


• First we set the parameters at very close and study the influence oi j3, f3 < 0, 
other parameters are fixed. Figlil illustrates the dependence of Hi{t) n f3. 

• Then we fix the parameter /3 (/3 = —2) and change an, the volatility of the 
first asset, other parameters are as in Figj^ Figj^the dependence of Hi{t) 
on (Til. Since 


H 




^1 

«'2Cr?i -t 4-3 ’ 


where 'I'l, 4^2,4'3 do not depend of an, if an increases, then the limit value 
Hi{t) becomes smaller. Thus, the asset with a small volotility is preferable. 
• Then we fixe (3 and an and study the influence of Ai. First of all we group 
the expression with respect to Ai and get 


TTOO _ 

■«1(3) - 


(q 2 - ^ 3 )^^! 
4>i 


(J>2 


where 4)i,<I>2 do not depend of Ai. Thus, the parameter Ai influences on 
the strategy for small t, whereas for large t this dependence is very weak 
provided at are close (see Fig§. 

The influence of an increasing of the risk parameter 7 analogous to an in¬ 
creasing of f3 by modulus. 


Let us summarize the influence of parameters on the character of the optimal 
strategy, analogous in the case of two and three assets: 

(1) an increasing of parameters /3 and 7 (by modulus) results a quicker attain¬ 
ment of the limit value as t —>■ (X); 

(2) for small time a decreasing of volatility of i-th asset (the values of aik) results 
an increasing of the proportion of this asset in the portfolio; 

(3) despite the fact that the trend Ai does not influence on the limit behavior as 
t ^ 00 , for small time the influence of this parameter if significant (increasing 
of Ai results the increasing of proportion of the corresponding asset). 


Remark 6.2. Let us note that for our strategy any moment of time can be taken 
as the initial one. Therefore, we can reasonably find the moment of time T for 
actualization of parameters of the model, i.e. for setting the time to zero. For every 
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set of parameters of the model there exists its own ’’infinity”, that is the time of 
achievement of the asymptotical value. For real data this time has an order of several 
years. It is natural to take this time as a time T for actualization. As it follows 
from our considerations, if the risk parameter 7 increases, then T becomes smaller, 
therefore, we have to actualize the model more frequently. 

Table 1 . Values of parameters for the case of two assets 



Fig 4 

Fig 5 

Fig 6 


-l;-5 

-1 

-1 

ai 

0.5 

0.5 

0.5 

0.2 

0.1 

0.1 

0.1 

7 

1 

1 

1 

Ai 

0.1 

0.1 

0.1;0.5 

^2 

0.1 

0.1 

0.1 

B 

1 

1 

1 

CTii 

0.1 

0.1;0.5 

0.1 

021 

0.1 

0.1 

0.1 

012 

0.1 

0.1 

0.1 

022 

0.1 

0.1 

0.1 

Ol3 

0.1 

0.1 

0.1 

O 23 

0.1 

0.1 

0.1 

Ai 

0.1 

0.1 

0.1 

A 2 

0.1 

0.1 

0.1 

A 3 

0.1 

0.1 

0.1 


Table 2. Values of parameters for the case of three assets 



Fig 7 

Fig 8 

Fig 9 

P 

-0.9;-2;-5 

-2 

-2 

Oi 

0.13 

0.13 

0.13 

02 

0.12 

0.12 

0.12 

Oz 

0.11 

0.11 

0.11 

7 

1 

1 

1 

Ai 

0.1 

0.1 

0.5;2;5 

A 2 

0.1 

0.1 

0.1 

A 3 

0.1 

0.1 

0.1 

B 

1 

1 

1 

<711 

0.1 
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7. A NONLINEAR INTEREST RATE (tHE COX-InGERSOLL-ROSS MODEL) 

7.1. An auxiliary problem: the conditional expectation and variance. The 

strict theory by Bielecki and Pliska is restricted to the case of the factor with a 
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Figure 4. Influence of parameter /3 on the strategy H = 

1. Hi, for P = -1; 2. H 2 , for /3 = -1; 3. Hi, for /3 = -5; 4. 
H 2 , for P = -5. 



Figure 5. Influence of parameter an on the strategy H — 
{Hi,H 2 ): 1. Hi, for an = 0.1; 2. H 2 , for an = 0.1; 3. Hi, for 
an — 0.5; 4. H 2 , for an — 0.5. 



Figure 6. Influence of parameter Ai on the strategy H = 
[Hi,H 2 ): 1. Hi, for Ai = 0.1; 2. H 2 , for Ai = 0.1; 3. Hi, 

for Ai = 0.5; 4. H 2 , for Ai = 0.5. 
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Figure 7. Influence of parameter (3 on the strategy 
(77i,7J2,i^3): 1. /3 = -0.9; 2. /? =-2; 3. /3 =-5. 





Figure 8 . Influence of parameter an on the strategy {Hi, H 2 , H 3 ): 
1. fJii = 0.3; 2. fjii = 1.5; 3. an = 3. 




Figure 9. Influence of parameter Ai on the strategy {Hi, H 2 , H 3 ): 
1. Ai = 5; 2. Ai = 2; 3. Ai = 0.5. 


constant volatility. The reason is that for this model the Hamilton-Jacobi-Bellman 
is reduced to a very special parabolic second order PDF where a sum of the order 
of derivatives and the order of polynomial in the coefficients at these derivatives is 
equal to two. In this section we consider another model of the interest rate where 
the volatility is proportional to the square root of the rate itself. The solution to 
problem can be found in this case for a special initial distribution of the interest 
rate. 


We consider a particular case of system (4.11: 

dF = {A + aR)dt + adWi, 
dR={B + j3R)dt + \VRdW 2 . 


(7.1) 

(7.2) 


Here H > 0, /3 < 0, cr > 0, A > 0, H, a are constants. 

The first equation describes the return F of asset with the trend that linearly 
depends on the interest rate R, that obeys the Cox-Ingersoll-Ross model [TU]. The 
inequality —2/3B > implies a positivity of the random process describing the 
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interest rate m- The interest rate of the form \7.2\ as a factor was considered, in 
particular, in the work [5] , a step toward to the finding the optimal strategy in the 
sense of Bielecki-Pliska [1]. Nevertheless the authors obtain only partial results. 

Let us assume that initially the interest rate R is distributed uniformly on the 
interval (0, L), L > 0. 

Remark 7.1. For the Vasicek model we obtain explicit formulae for initial Gaussian 
distribution including limiting cases. Nevertheless, for the Cox-Ingersoll-Ross we 
can get an explicit formula only for uniform initial distribution. 

The Fokker-Planck equation for the join distribution of random values F and R 
given by system ( [7T| ), is 

mt. f, r) ^ ^ + ,5P„,/,,)+ 


dt 

+ {B + /3r - A^) 


df 

dP{tJ,r) 1 ,d^P{tJ,r) 1 2 d^PitJ,r) 


dr 


2 ^ 


dp 


— -A^r= 
2 


= 0 , 


(7.3) 

subject to initial conditions 

P{0,f,r)=Hf-fo)X(o,L)(.r). (7.4) 

Rigourously speaking, to define a probability density function P (0, /, r) we have 
to divide the expression (7.4) by L. Nevertheless, as follows from the linearity of 


equation (7.3) and definitions (4.17), (4.18) this multiplier does no influence on the 


result of computations. 

The Fourier transform with respect to (/, r), the function P(t,/r,^), obeys the 
equation 

+ ( A + Bp] P {t: 7^1 0 = 0, 


dt 


(7.5) 


with initial conditions 
1 


m = ^ ^ ^ P Sp) as L ^ oo. (7.6) 

27r 27r 


Equation (7.5) has the first order and can be integrated. The solution to the 
problem (7.5),(|7.6|) in the limit case L —>■ oo can be found by a standard way: 


^ 2 i/oM a 2 + 2 t A^l i + 2 t B l3-\-t 

P{t,fJ-A) = e~ 2 X 2 S{s{t,^i,p)x 

X'^{2ia + 2iP ^ + X^p) cosh ( ^ JifL _j_ ^ arctan( , ^ =)^ 

y Y'2 zq;/x A2+/32 j 


2ia^X‘^ + 




where 


si.t,n,p = - 


(^p — p2ia~fiXFppp+2 a fije ^_2 a/i 


p2ia fi A 2 +/ 92 _|_j^ 2 ^_^_|_ ^.^2ia^A2+/52—jA2^+/3^ g-£v77i,yyA2+^ 


We apply (4.17) and (4.18): 


Pit,O,p = 0{t,pS{sit,O,p), 

d^P{t, 0, p = pit, p Spit, 0, p) + 'Pit, p Spsp, 0, p). 


(7.7) 

(7.8) 


dlPit,0,P = qpt, P Spit, 0,p) + q 2 it,p d' isit,0,p) + q3it,p6Psit,0,p), (7.9) 
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where 






e{t,0 = 


+ (2/3 — ’ 

-/A^(*^^A^ — 2/3^) cosh + i arctan(^-^t^)^ 


Li — 


Ht, 0 = 0 {Li - i{At + fo )), '0(t, 0 = 0{t, OL 2 , 

Ba ((4/3A4^2 _ (4^2 _ 2A4^2 _ 6/^A2^)) 


X“ 


/3A2^(A2^ + 2i/3)2 

sinh(f+zarctan(^^)) 2ag(/A"e(A^e" - 5^") - A/SX^e) + 2/3^ 

cosh(^+/arctan(A!^)) ^ + 2iPY 

-a^2;^4g2/3t 2((A2^ + 2i/3)X^^t - 2iX^^ + 2/3)aPe^* + 4ial3X^^ + aX^^ - Aa/3^ 

//3A2^(e/3* - 1 ) + 2/32 ■ 

We do not write here the explicit values of qi{t,^),i = 1,2,3, since they are very 
long. 


J 2 — 


We substitute (7.7),(7.8),(7.9) in (4.17) and (4.18) and after cumbersome but 


(7.10) 


standard computations get 

- / a A (1 - e^‘)ar (2 + A^)ae^^* 

(l+A^)at (2+A^)a\ 

P - 1-^2 ) ^ 

^ ( 2e^P* (2/3t + 3)e2'5* 2e^* lA 2 x 2 

=ta + +- ^3 -^ ^ 

/ 4 ., 3(4/3f + 5)e3/5* ^ ^ ( 2 / 32 ^ 2 -2/3t-5)e^* A a2A4 

+ ( ^ ^ ^^ 


_ 2B(^« + 3)P« + - 4fie« - 


2 x 2 


a^A 


(7.11) 


7.2. Example of portfolio consisting of two assets. We consider a simple as¬ 
set allocation example, featuring an interest rate which affects a stock index and 
also serves as a second investment opportunity, illustrates how factors which are 
commonly used for forecasting returns can be explicitly incorporated in a portfolio 
optimization model. This example was systematically considered in the works by 
T.Bielecki and S.Pliska (see i, IS], i). The dynamics of the security prices is 

dS'i(t) 


Slit) 


= {Ai -I- aiRit))dt + aidWi{t), 


dS2it) 

S2it) 


= Rit)dt, 


where R{t) is the Cox-Ingersoll-Ross interest rate (7.2). 

dV ' dS dS 

The capital of portfolio obeys the equation — = h— —I- (1 — h)-—, where the 

V bl 02 

scalar valued function h is interpreted as the proportion of capital invested in the 
risky asset, leaving the proportion 1 — h invested in the bank account. 
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Remark 7.2. We can consider a portfolio consisting of any number of assets, as it 
was done for the case of the Vasicek-type interest rate. 


Let us denote InV = F. Then 


dF = ^Aih — + (^{h — 1) «! + dt + aihdWi. (7. 


12 ) 


For the system (|7.12|), (|7.2|) we apply the formulae for conditional mathematical 

(7.13) 


expectation and variance (7.10), (7.11) after substitution 


1 


A = Aih—^alh^, a = (ai — l)/i + 1 , a = aih. 


We consider again the functional ( |5.2[ ): 

Q^{t, r; h) = f{t, r; h) - 'yv{t, r; h), 

where 7 is the risk aversion coefficient, and find the optimal strategy in the sense of 
Definition 15.11 

According to (7.10), ( 7.11| ) and (7.13) we get 

Qj{t, r; h) = K 2 h? + Kih + Kq, 

where Ki and smooth functions of t,r and coefficients Ai, ai, cri, il,/3, A, 7 . These 
functions can be expressed through elementary functions, nevertheless, these expres¬ 
sions are cumbersome and we do not write them. Since Q.y{t, r; K) is quadratic with 
respect h, and 

ta^h? 

AI2 ^ ^)l(a—ai —1, a—cr\) ^ O5 

then Q.y(t,r]h) has a unique point of maximum (analogous to the linear case, see 
Sec|^, the respective optimal strategy in the sense of definition [5^ is the following: 

-Ki 


= 


2 K 2 


(1 - -k -f Mie^* -k Mq) + AifdH 

(1 - ai)2(M4e4/3‘ -p M3e3/3t + N2e'^P* + + Nq) + {2j + l)all3H' 


where 


M 4 = jX^{8X^ + 55), M 3 = -jX^ ((3A2 -f B)Al3t + 15X^ + 4r/3 -k 125) , 

M 2 = 2-fX‘^ph'^ + (2A2 -f /3r)4jPX^t + 127A^ -k (55 -h 2/3r)3'yX'^ - (X^ + B)/3^, 
Ml = 2-iX'^(3H'^ + (^2 _ 2-fX^)l3X‘^t - 5jX^ + {(3^ - - 8'fB)X^ + /3V -h P^B, 

Mo = {P‘^-2'yX^)/3Bt+2jX^Pr-P^r, Ni = -fX^{2(3'^X'^t^-2(3X^t-5X^-4/3r-8B), 
N 2 = ')X^{2X^0^t^ + {2X^ -f /3r)4/3t -k 12X^ F 6l3r + 155), Nq = 2jl3X‘^{-Bt + r). 
We get 


lim 5.y = lim H.y{t) = 


_ (ai - l)(/3^5 - 2-fX^B) - Ai/g^ 

t->-oo t-loo ^ (ai — 1)227A^5 -|- (27 -|- l)af/3^ 

As was shown in Sec|^ in the case of a linear interest rate model 


lim H.y(t) = 


-1 

«! — 1 ’ 


ai 1 - 


(7.14) 


This expression can be obtained by a limit pass in (7.14) as 7 —)■ 00 and cti —> 0. 
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8. Comparing the optimal strategies for linear and nonlinear models 

OF THE INTEREST RATE 

We set the interest rate r = 0.05, the initial capital of portfolio /o = 0.08, the risk 
aversion coefficient 7 = 0.1, the parameters B = 0.05, /3 = —1, A = 0.04, Ai = 0.15, 
ai = — 1 , tJi = 0.2 are taken from 

Fig. [^illustrates the corresponding optimal strategies of investment for the case 
of linear (solid line) and nonlinear (dashed line) interest rate for their uniform initial 
distribution. We see that these strategies are very different. To approach the optimal 
strategy for the nonlinear case to the strategy for the linear case we should choose 
A larger risk sensitive parameter 7 . 



Figure 10. Comparing the optimal strategies for different models 
of the interest rate: 1. Vasicek-type interest rate 2. Cox-Ingersoll- 
Ross interest rate. 


As we have seen, in the case of the Vasicek-type interest rate the asset less de¬ 
pendent on the factor is preferable for the investment for a large time. As follows 
from different combination of parameters, for the Cox-Ingersoll-Ross interest rate 
the properties of the factor are taken into account more effectively. 


Remark 8.1. If we assume additionally that the variance of the return of the risky 
asset satisfying ( |7.1 1 is proportional to the interest rate, we fall in the situation of 
the Heston model |iTJ|. one of the most popular models of stochastic volatility. In 
|30j we analyze the value of mean dispersion, and formula (4.171 turns useful there, 
too. 


Remark 8.2. Papers m and [IH] extend the work [5], they deal with the problem 
of long-run optimal investment in the frame if the Cox-Ingersoll-Ross model. 

Remark 8.3. As follows from [2], the behavior of strategy of investment should 
depend on the speed of decay at infinity the initial distribution of the interest rate. 

Remark 8.4. Equations ( fni ), refer to the so called “affine” model therefore 
the respective Fokker-Planck equation can be solved explicitly. 
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